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The transport properties of molten LiF-YF3 mixtures have been studied by pulsed field gradient nuclear 
magnetic resonance spectroscopy, potentiometric experiments, and molecular dynamics simulations. The cal- 
culated diffusion coefficients and electric conductivities compare very well with the measurements accross 
a wide composition range. We then extract static (radial distribution functions, coordination numbers dis- 
tributions) and dynamic (cage correlation functions) quantities from the simulations. Then, we discuss the 
interplay between the microscopic structure of the molten salts and their dynamic properties. It is often 
considered that variations in the diffusion coefficient of the anions are mainly driven by the evolution of 
its coordination with the metallic ion (Y 3+ here). We compare this system with fluorozirconate melts and 
demonstrate that the coordination number is a poor indicator of the evolution of the diffusion coefficient. In- 
stead, we propose to use the ionic bonds lifetime. We show that the weak Y-F ionic bonds in LiF-YF3 do not 
induce the expected tendency of the fluoride diffusion coefficient to converge toward the one of yttrium cation 
when the content in YF3 increases. Implications on the validity of Nernst-Einstein relation for estimating the 
electrical conductivity are discussed. 
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I. INTRODUCTION 

Molten salts are liquid mixtures of ionic species at 
high temperatures. They receive much attention due to 
their diverse but promising applications that range from 
metal production*—, the pyrochemical treatment of nu- 
clear wasted—, primary or secondary coolants in several 
concepts of generation IV fission 8 and fusion 9 nuclear re- 
actors, and electricity storage devices 1 ^. 

In the context of nuclear applications, the nuclear reac- 
tions produce rare earth elements in the core. From the 
neutronic point of view, these have a poisoning effect and 
must therefore be eliminated. One of the advantages of 
most of the concepts of reactors based on molten salts, in- 
cluding the promising thorium-fueled molten salt fast re- 
actor (MSFR) is that their fuel is dissolved in the primary 
coolant (which is the molten salt) 11 ' 12 . Such liquid state 
fuel allows for the elimination and reprocessing of sev- 
eral fission products and rare earth elements during the 
reactor operation, via chemical or electrochemical meth- 
od o 13 ' 14 . One drawback is the complex handling of the 
molten salt fuel inside the (recycling) circuit: important 
care must be taken for controlling the physico-chemical 
condition a 15 ' 16 , and additional efforts have to be made on 
the development of adapted structure material a 8 ' 17 ' 18 . 



a ) Electronic mail: maximilien.levesque@gmail.com 
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The knowledge on molten fluoride salts chemistry has 
long remained fragmented and limited to a few experi- 
mental techniques (mainly electrochemical ones) due to 
experimental hindrances. One has to deal with high tem- 
peratures (from 500 K to 1800 K) and to control the 
corrosive properties by setting up original experiments. 
Many efforts have been made towards this direction and 
several experimental techniques have been adapted to the 
study of molten salts. The development of new cells has 
unlocked the study of molten fluorides by EXAFS spec- 
troscopy^, which, coupled to other spectroscopic meth- 
ods such as NM Bi 20 ' 21 or Rama n 22 ' 23 , is able to provide a 
view of the local structure in these melts. As for the dy- 
namic properties, similar improvements have been made 
and the self-diffusion coefficients can be measured up to 
1500 K by pulsed field gradient nuclear magnetic reso- 
nance (PFG-NMR) 24,25 . Such measurements are impor- 
tant since the diffusion coefficients provide an informa- 
tion for the individual dynamics of the ions, complement- 
ing the viscosity or electrical conductivity measurements 
which inform on the collective dynamics of the macro- 
scopic system. 

During the past decade, these experimental efforts 
have been accompanied by the development of simula- 
tion tools adapted to the study of molten fluorides and 
chlorides, namely density functional theory 2 ^, thermody- 
namic modeling 2 ^— and molecular dynamics (MD) 30 ' 31 . 
Among those, the latter has proven to be very efficient 
for the study of the structure and dynamics of the sys- 
tem at the microscopic scale, allowing for the interpre- 
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tation of experimental results^ — and the prediction of 
quantities that hitherto remain unknown^. The simu- 
lations are based on the use of a polarizable ion model, 
which has strong physical grounds^— and can efficiently 
be parameterized from first-principles^—, i.e., without 
having empirically recourse to any experimental data. 

In this work, we combine state-of-the-art experimental 
and simulation tools to study a series of LiF-YF 3 molten 
salts of varying compositions. After describing shortly 
the methodological aspects, we discuss measured and cal- 
culated self-diffusion coefficients and electrical conduc- 
tivity. We show that MD simulations are in quantitative 
agreement with experiments, which allows for further use 
of the simulations to discuss the interplay between the 
microscopic structure of the electrolytes and their dy- 
namic properties. 

II. METHODS 

A. Pulsed field gradient nuclear magnetic resonance 

Mixtures of LiF and YF 3 salts (purity 99.99 %) were 
prepared in a glove box under argon in order to avoid 
H2O and O2 contamination of the samples. We use boron 
nitride (BN) crucibles to confine the salts, and to prevent 
them from any chemical reactions. The amount of salt in 
each crucible is ca. 50 mg. The BN crucible is heated by 
a symmetrical irradiation with two CO2 lasers. The laser 
power is increased progressively until the target temper- 
ature is reached. 

The high temperature PFG NMR spectra were 
recorded using a Bruker Avance WB 400 MHz operating 
at 9.40 T. An argon stream prevents oxidation of the BN 
at high temperature. The 10 mm liquid NMR probe is 
specially designed by Bruker and adapted in CEMHTI to 
work up to 1500 K. It is equipped with a gradient coil pro- 
viding 5.5 G-cm -1 A -1 that is combined with a gradient 
amplifier of 10 A. We use a NMR pulse sequence combin- 
ing bipolar gradient pulses and stimulated echo 4 ^. It is 
repeated with 8 gradients of linearly increasing strength. 

Measurements were performed using the following 
NMR parameters: durations of radiofrequency magnetic 
field application for ir/2 pulses p90 = 13 /lis ( 19 F) and 
p90 = 9 /js ( 7 Li), gradient strength g was varied from 
to 50 G/cm, gradient application time S between 1 and 
5 ms. For 7 Li measurements, we used a pre-saturation 
cycle before the diffusion sequence because of the long 
relaxation time T± of this nucleus (25 s in pure molten 
LiF). 

Finally, the self-diffusion coefficients were obtained by 
nonlinear least-square fitting of the echo attenuation 4 ^. 

B. Electrical conductivity potentiometric measurements 

The experimental setup used for measuring the electri- 
cal conductivity is described in details in reference 4 ^; the 



conductivity cell configuration was derived from the one 
proposed by Hives and Thonstad™. The conductivity 
was measured for LiF-YF3 mixtures with compositions 
ranging from xyp 3 — 0.00 to 0.18, at several tempera- 
tures. Here we report the results obtained for the 1130 K 
isotherm. All the details are given in the Supplementary 
Material Document No.l— . 



C. Molecular dynamics simulations 

We use molecular dynamics (MD) to compute static 
and dynamic properties of the molten salt at the atomic 
scale. The MD simulations used therein rely on inter- 
atomic interaction potentials parameterized from first 
principles. The parameterization procedure is described 
elsewhere^. They are of ab initio accuracy and have 
been successful in reproducing accurately the electrical 
conductivity, the heat capacity, the viscosity, etc, of var- 
ious ionic materials (mainly molten fluorides and solid 
oxides 45 ' 46 ). In the case of molten fluorides, the inter- 
action potential is best described as the sum of four dif- 
ferent components: charge-charge, dispersion, overlap re- 
pulsion, and polarization. The charge-charge term is nat- 
urally 

V™(r<,)=£^, (1) 

i.J>i 

where q+ is the formal charge on ion i and r\j is the dis- 
tance between ions i and j. The dispersion component 
captures dipolc-dipole and dipole-quadrupole terms. It 
reads 




where Cq (Cg J ) is the dipole-dipole (dipole-quadrupole) 
dispersion coefficient, and /g J 8 are Tang-Toennies disper- 
sion damping functions^ 7 -. They describe the short range 
correction to the asymptotic multipole expansion of dis- 
persion 48 . They are expressed as 

/« (r l0 ) = 1 - 4 exp (-b% rij ) £ (3) 

fe=0 

where Vfi is the distance at which the correction starts 
to be taken into account, and cjf = 1, ensuring that 
the dispersion term is cancelled for short distances. The 
repulsion overlap component is given by 

V tep (r ii )= ^A«exp (4) 

i,j>i 

where A lJ and a lJ are also parameters deduced from first 
principles. Finally, the interatomic potential includes 



charge-dipole and dipole-dipole polarizations: 
V p ° l (n 3 ) = E (;/,/< ...„/!' (r«) - r/,,/,.,,/T (r«) *) T« 
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Here, and ^ are the charge-dipole and dipole- 

dipole interaction tensors, and ccj is the polarization of 
ion i. Tang-Toennies functions are again included to ac- 
count for the short-range effects. The induced dipoles fii 
are three-dimensional vectors. At each time step, they 
are relaxed adiabatically and self-consistently in order to 
minimize V p . The inter atomic potential parameters 
are given in reference 3 ^. 

We studied 11 concentrations ranging from 
xyf 3 = 0.10 to 0.60 by steps of 0.05. These sys- 
tems have been simulated along two isotherms. The 
first one at 1130 K enables comparison with experi- 
mental data. The other one at 1200 K for comparison 
with previous work on fluorozirconate a 34 ' 49 . For both 
temperatures, (i) the equilibrium volume was deter- 
mined from a simulation in the NPT ensemble with a 
target pressure of 1 atm, using the method described 
by Martyna et al£Q; (ii) production runs were then 
performed in the NVT ensemble^ with a time step of 
1.0 fs. The total simulation time was of 4 ns for the 
production runs for each system. The relaxation time 
of the thermostat was set to 20 ps. The simulation cell 
parameters and the number of atom of each species are 
summarized for each composition and temperature in 
the Supplementary Information^ 4 -. The cut-off distance 
for the real space part of the Ewald sum£2 and the short 
range potential was set to L/2 where L is the length of 
the cubic simulation cell. The value of the convergence 
parameter in the Ewald sum was set equal to 5.6/ L 
and 8 3 fc-vectors were used for the reciprocal space part 
calculation of the forces and potential energy. 



III. TRANSPORT PROPERTIES 

The self-diffusion coefficients of fluoride, lithium and 
yttrium ions along the 1130 K isotherm are presented in 
Fig. [TJ as a function of the salt concentration. Both ex- 
perimental and molecular dynamics results are presented 
for F~ and Li + . For Y 3+ , the self-diffusion coefficient can 
only be measured from molecular simulations, due to the 
NMR characteristics of 89 Y (low gyromagnetic ratio and 
long relaxation time up to 45 min). The self-diffusion 
coefficient D a of species a is extracted from molecular 
dynamics simulations using Einstein's relation, textiti.e., 
from the long-time slope of mean squared displacement: 



D(a) = l lim d t (\ n (t)-r t (0)| 2 ), 



(6) 



where r; (t) is the position of atom i of type a at time 
t. The brackets denote an ensemble average over all tra- 
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FIG. 1. Self-diffusion coefficients of lithium (black), fluoride 
(red) and yttrium (blue) ions at 1130 ± 10 K as a function of 
the salt concentration. Full (open) symbols indicate experi- 
mental (Molecular Dynamics) data. Dashed and dotted lines 
are guides for the eyes. 



jectories and atoms i S a. The self-diffusion of Li + is 
about 60 % higher than that of F~, whatever the con- 
centration. Excellent agreement is found between exper- 
iments and simulations. We can therefore be confident 
in the MD simulation results for Y 3+ . The latter value 
is about one third of the diffusion coefficient of Li + over 
the whole range of compositions. All self-diffusion co- 
efficients decrease at the same rate with the increase in 
yttrium fluoride concentration. 

The electrical conductivities measured at a temper- 
ature of 1130 K are presented in Fig. [5] for different 
salt compositions, with previous experimental results ob- 
tained for larger YF3 concentrations^ 3 -. In molecular dy- 
namics simulations, this quantity is calculated from 



V 



lim 



N, , 



^q a A a (t) 



(7) 



where j3 — 1/kgT with fc^ the Boltzmann constant, e is 
the elementary charge, V the volume of the simulation 
cell, and A a (t) = J2iea & r i (t) is the net displacement of 
all the ions of species a and charge q a in time t. There 
are N a — 3 species here. The error bars associated to 
our calculations of the electrical conductivity from Eq. [7| 
is estimated to 10 %. Such large error bars are due to the 
collective nature of \- It is n °t averaged over both the 
time steps and all the atoms of a given species, as self- 
diffusion coefficients here. Therefore, one often approxi- 
mates calculation of Eq. [7] by using the Nernst- Einstein 
relation, which relates x to the individual diffusion coef- 
ficients of the various species. Surprisingly, we find the 
electrical conductivity calculated from Nernst-Einstein to 
be overestimated by 10 to 20 % only (larger differences 
are usually observed for such salts). MD results from 
Eq. [7] are plotted in Fig. [5J Good agreement is found 
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FIG. 2. Electrical conductivity % of molten LiF-YF3 mix- 
tures as a function of the salt concentration at 1130 K. (black 
circles) experimental data from this study; (white circles) ex- 
perimental data from reference—; (red squares) molecular dy- 
namics results. 

again between experiments and simulations. 

The electrical conductivity falls from s» 9 r2 _1 cm" 1 in 
pure LiF to « 3 fi cm -1 , when the molar fraction of YF3 
rises from to 0.5. This decrease in \ despite the addition 
of the highly charged Y 3+ ions is not surprising in view of 
the variations observed in other mixture o 30 ' 32 ' 53 . It cor- 
relates with the reduced diffusion of all species with iyf 3 
observed above; in the following sections, we will analyze 
these variations in light of the microscopic structure of 
the melt. 



IV. MICROSCOPIC STRUCTURE 

In the previous section, we have shown that our molec- 
ular dynamics simulations were able to predict the trans- 
port properties of molten LiF-YF3 with a very high ac- 
curacy. An advantage of such simulations is that we can 
extract precisely atomic scale informations on the local 
structure of such melts. We will mainly focus on the 
study of the solvation shell of yttrium ions. 

In Fig. [3J we present the partial radial distribution 
functions (RDF) gr y (r) for F--F", F~-Y 3 +, Y 3 +-Y 3 + 
and Li + -F _ pairs at 1130 K. Situations corresponding to 
low and high contents in YF3 (xyf 3 = 0.15 and 0.60) 
are presented. Firstly, we observe that the Y-F ionic 
bond length, taken at the position of the first maximum 
in g(r), is of «2.24 A for all the compositions. This is 
significantly larger than the Zr-F distance («2.02 A) in 
fluorozirconate melts^ 4 ^ while both cations have rather 
similar short-range repulsion potentials^. This is due 
to the larger charge carried by the Zr 4+ cations which 
therefore tend to be more attractive to fluorides with 
opposite charges. Secondly, we observe that the first peak 
in this Y-F RDF is twice less intense and sharper than 




r(A) 

FIG. 3. Radial distribution functions at 1130 K for iyf 3 = 
0.15 (black line) and 0.60 (red dashed line). 



its Zr-F counterpart. Also, the following minimum is less 
marked in YF3. These quantities seem to show that the 
first solvation shell of Y 3+ is more weakly bound than 
the one of Zr 4+ in molten fluorides. In the case of the 
Li + -F~ pair, the RDF has a shape very similar to the 
case of other molten fluoride mixture^. 

The extraction of precise coordination numbers is a dif- 
ficult experimental task, for example by extended X-ray 
absorption fine structure (EXAFS) spectroscopy, that 
leads to errors of the order of 20 % for liquid state sys- 
tems. Here the difficulty is increased by the high tem- 
peratures of the systems. From NMR, it is possible to 
extract an average coordination number by comparing 
the chemical shift of the probed species to typical values 
obtained for this quantity in well-defined solid state com- 
pounds. This was performed in LiF-YF3 mixture o 54 ' 55 , 
which led to the conclusion that the yttrium ions are 
most probably 7-fold or 8-fold coordinated. A smaller 
coordination number (CN) of 6 was extracted from Ra- 
man spectroscopy experiments by Dracopoulos et al. in 
KF-YF3 mixture o 56 ' 57 . Such a reduction in the CN when 
switching from LiF to KF-based systems and a stabi- 
lization of octahedral coordination shell in the latter is 
consistent with our recent results in molten fluorozir- 
conates^ 4 -. 

Because of the intrinsic atomic-scale of our MD sim- 
ulations, we can extract CNs with estimated error less 
than 2 %. The evolution of the CN of yttrium with its 
molar fraction in LiF-YF3 is presented in Fig. 01 We 
find the predominant species to be 7-fold coordinated yt- 
trium ions for xyf 3 below 30 % and 8-fold coordinated 
ions above. From NMR experiments, one extracts the 
average CN over time. From Fig. [H the average CN for 
an yttrium ion at the eutectic composition (xyf 3 = 0.2) 
is « 8.3 in our simulations, which compares well with the 
one estimated from NMR experiment o 54 ' 55 . We observe 
almost linear evolutions for all CNs. The proportions of 
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FIG. 4. Evolution of the calculated coordination numbers 
with the composition of the melt. CNi stands for the pro- 
portion of i-fold coordinated ytrium ions. No noticeable dif- 
ferences were observed between the two temperatures studied 
here (1130 K and 1200 K). Reported values have a typical 
error bar of 1 %. 



CNs 6 and 7 decrease with a;YF 31 and are compensated 
by increases in CNs 8 and 9. Contrarily to the case of 
LiF-ZrF4 mixtures, no anomalous behaviour is observed 
around a given specific composition 49 . 

As for the like-like RDFs which are also shown in 
Fig. [3j the F-F one does not show any strong evolution 
with the composition of the melt. The Y-Y one is char- 
acterized by a first peak, which is located at a distance 
of «3.9 A which is smaller than twice the Y-F distance 
and therefore corresponds to pairs of cations which are 
sharing a common fluoride anion. The intensity of this 
first peak increases with the YF3 concentration, and it 
is larger than the one of the corresponding Zr-Zr RDF is 
fiuorozirconate a 34 ' 42 ' 49 . Such a structural feature, which 
denotes some medium-range, network-like ordering, has 
been shown to impact importantly the transport proper- 
ties in molten fluorides. 



V. DISCUSSION: INTERPLAY BETWEEN STRUCTURE 
AND DYNAMICS 

In molten salts, the coordination number is often seen 
as the most important structural criterion. It is heavily 
used for the interpretation of various properties through 
the perspective of the Lux-Flood acidit y 58 i 59 . In a molten 
fluoride melt, this acidity's definition lead to an acid be- 
ing a fluoride ion acceptor and a base being a fluoride 
ion donor—. In this context, based on their coordina- 
tion numbers, which range between 6 and 9, Y 3+ and 
Zr 4+ should both be considered as strongly acidic species. 
Nevertheless, very distinct behaviours are observed in 
these melts concerning the transport properties. 

On the one hand, in MF-ZrF4 melts (with M + be- 



ing Li + , Na + , K + or even a mixture of Li + and Na + ), 
the Df /Dm ratio starts from a value very close to unity 
in the alkali fluoride melt, then strongly decreases upon 
addition of ZrF4. In parallel, the Dp/Dzr ratio also 
decreases, taking values close to unity in ZrF4 concen- 
trated melt o 34 ' 61 (note that a similar behavior is observed 
in LiF-BeF 2 mixtures 3 ^^). Here, in L1F-YF3 both ra- 
tios hardly vary, showing that the diffusion of fluoride is 
barely impacted by the presence of yttrium ions, even at 
high YF3 concentrations. 

On the other hand, we have already mentioned that the 
Nernst-Einstein approximation overestimates the electri- 
cal conductivity of L1F-YF3 molten salts by 10 to 20 %, 
independently of the composition. Such a relatively good 
agreement is usually observed in molten alkali halides, 
textiti.e., in the absence of a complexing species. In 
LiF-NaF-ZrF4 mixtures, we have shown that deviation 
from the Nernst-Einstein approximation increases with 
xzrF 4 , the content of complexing species. For the highest 
concentration studied (xz r F 4 =0.29) the conductivity is 
overestimated by as much as a factor of 2 if correlations 
between the ions are neglected. This again shows that 
association effects between Y 3+ and F~ ions are weak in 
L1F-YF3 melts, despite the high coordination numbers. 
These points lead to the picture of a loosely bound net- 
work structure. The association behaviour of L4F-YF3 
is thus qualitatively different from that of LiF-ZrF4, al- 
though their coordination number is high. 

In order to investigate further the associative be- 
haviour in LiF-YF 3 , we analyze the rate of change of 
the instantaneous coordination number, which can be in- 
terpreted as the solvation shell lifetime. This is done by 
performing a cage correlation function analysis. The first 
minimum of the Y-F partial radial distribution function 
is used as a simple geometrical criterion to determine the 
limit of the first coordination shell of a given Y 3+ . The 
number of counter-ions F~ that have left the shell be- 
tween times t and t + St is noted n° ut (t,t + 6t). The 
rate at which a single counter-ion (F - ) leaves the Y 3+ 
ions solvation shell (the cage) is given by the cage-out 
correlation function 62 , which is defined as: 



cr (*) = (e(i-nr t (o,t))) 



(8) 



where O is the Heaviside step function. 

The cage-out fonctions for a =Y 3+ and Zr 4+ in LiF- 
YF3 and LiF-ZrF4 are presented in Fig. [5] for the same 
molar fraction of 0.35 at 1200 K. Note that the changes 
with composition are very small, as shown in the supple- 
mentary materials. The exponential decay of the cage- 
out function is much faster for the Y 3+ -based salt. The 
corresponding characteristic time is of 0.7 ps, which is 
much lower than the ones obtained in LiF-ZrF4 (3.1 ps), 
NaF-ZrF 4 (15.7 ps) and KF-ZrF 4 (76.4 ps). The very 
short characteristic time measured in the case of LiF-YF3 
melts is at the origin of the particular transport proper- 
ties obtained in this melt. It shows that even if Y 3+ and 
Zr 4+ have rather similar (and high) coordination num- 
bers, the ionic bond is much weaker for the former. 
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FIG. 5. Cage-out correlation functions for Y 3+ and Zr 4+ 
ions solvation shell in LiF-YF3 and LiF-ZrF4, respectively, 
at 1200 K and molar fraction iyf 3 or xzvF 4 = 0.35. 



Although the use of such cage correlation functions 
appears to be very useful for understanding the associa- 
tive properties of molten salts, this quantity can only 
be obtained from simulations. It would nevertheless be 
valuable to extract similar information from a quantity 
which is more readily accessible. To this end, we com- 
pare in Fig. [6] the effective potential between Zr-F and 
Y-F ions in the two melts, which was obtained from: 



pW??{r) = -hx 9ij (r)-2\n{r). 



(9) 



Indeed, Eq. [9] relates the partial radial distribution func- 
tions gij (r) discussed above to an effective interaction po- 
tential. The last term in Eq.|9] accounts for the rotational 
invariance of Wf^{r). Following the transition state the- 
or y 63 ' 64 , the rate at which the F~ ions escape from the 
Y 3+ solvation shell should depend strongly on the height 
of the corresponding barrier. We observe a smaller bar- 
rier (by k, 1 kgT) in the LiF-YF 3 melt, in agreement 
with the smaller characteristic time extracted from the 
cage-out correlation function. This clearly shows that 
the shape of the RDF, rather than the CN distribution, 
should be studied for determining the associating prop- 
erties in a molten salt. 
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have shown that despite being appealing, the coordina- 
tion number is not the right observable to predict the 
strength of the association between a metallic species and 
the anion. It does therefore not bring much insight on the 
transport properties of molten salts. Indeed, although 
Y 3+ and Zr 4+ ions have similar coordination numbers in 
fluorozirconate melts, their diffusion coefficients do not 
show similar behavior: In LiF-ZrF4 melts, when a;zrF 4 
increases, Dp decays much faster than Dn, eventually 
reaching £>z r . On the contrary, in L1F-YF3 melts, the 
three diffusion coefficients evolve in a similar way. 

The analysis of the effective potentials, which can read- 
ily be extracted from the radial distribution functions, 
and of the cage correlation functions, which measure the 
lifetime of the first solvation shell of the ions, provide a 
much more complete framework for interpreting the dy- 
namic properties. Here we show that the variation of the 
diffusion coefficients is linked to a short lifetime of the Y- 
F ionic bonds, which in turn arises from a rather shallow 
well in the effective potential. This points toward the po- 
tential use of the transition state theory for rationalizing 
the ionic association effects in molten salts. 



VI. CONCLUSION 
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